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UNDERSTANDING TRANSITION AND TURBULENCE 
THROUGH DIRECT SIMULATIONS 


P. R. Spalart &: J. J. Kim 

A. R. C., C. F. D. Branch, Turbulence Physics Section 


Direct simulations consist in solving the full Navier-Stokes equations, 
without any turbulence model, and describing all the detailed features of the 
flow. Usually the flows are three-dimensional and time-dependent and contain 
both coarse and fine structures, which makes the numerical task very chal- 
lenging in terms of both the algorithm and the computational effort. Most 
of the work until now has involved spectral methods, which are highly accu- 
rate but not very flexible in terms of geometry or complex equations. For 
that reason, future work will also rely on high-order finite-difference or other 
methods. 

Direct simulations complement experimental work, and both contribute 
to the theory and the empirical knowledge of turbulence. Once such a sim- 
ulation has been shown to be accurate the flow field is completely known, in 
three dimensions and time, including the pressure, the vorticity and any other 
quantity. On the other hand, most simulations to date solved the incompress- 
ible equations in rather simple geometries, and direct simulations will always 
be limited to moderate Reynolds numbers. Extensive simulations have been 
conducted in homogeneous turbulence, channel flows, boundary layers, and 
mixing layers. Much effort is devoted to addressing flows with compressibil- 
ity and chemical reactions, and to new geometries such as a backward-facing 
step. 
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UNDERSTANDING TRANSITION AND TURBULENCE 
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Lyapunov exponents of a turbulent channel flow. 
Growth of compressible mixing layers. 


Secondary instability in channel flow 
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Evolution of Modal Energy 


si 



,.ovot.o;.o;.oi.oi.o;.oi.oj,.oA oj, .qj oA.QA.qj.oi 


140 


0.0 600.0 1200.0 1800.0 2400.0 3000.0 3600.0 4200.0 4800.0 5400.0 6000.0 


Stability and relaminarization of the swept attachment-line flow 
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(i^, jR) space, where K is the suction parameter and R the Reynolds 
number based on strain rate and sweep velocity component. Suction 
has much less effect on relaminarization than on linear stability. 



Critical Re at Leading Edge 





Pressure terms in a Reynolds-stress turbulence model 
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• As a result, the widespread k — e and Reynolds-stress models incur 
large errors near walls. 



Pressure-strain term 
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Vtuv KK '21 . Pressure strain term, p'u', 2 , in the budget equation for u' 2 a 2 across the channel. O term 
computed from the channel data; . model. (equations (4ti) + (48) + (4!>)) ; . model, 

equation (4(i). . model, (equations (4H) + (411)). 


Coherent structures in shear flows 
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• When S* is given such values in a homogeneous shear flow, streaks 
appear, indicating that the presence of a wall is not essential to form 
them. 
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The dimension of strange attractors in turbulent shear flows 
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complexity of a fully developed turbulent now. 

• This result shows that shear flow turbulence cannot be considered 
to result from the interaction of a “few” degrees of freedom. 
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• S imu lations have: a) reproduced this effect; b) validated M c ; 
c) provided a physical explanation. 



